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Abstract. Strong final-state interactions create a pronounced cusp in r/' 
generalize the non-relativistic effective field theory framework developed for the extraction of tttt scattering 
lengths from K ^ Stv decays to this case. The cusp effect is predicted to have an effect of more than 8% 
on the decay spectrum below the 7r~^7T~ threshold. 
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1 Introduction 

In the last few years, the investigation of the cusp ef- 
fect in the decay 7r^7r^7r+ has become one of 
the most precise methods to extract S-wave pion-pion 
scattering lengths from experiment [1-6]. Very loosely 
speaking, the cusp in the invariant mass spectrum of the 



r- fol- 



TT^TT^ pair is generated by the decay 
lowed by charge-exchange rescattering tt+tt" 
plus the fact that the pion mass difference shifts the 
7r+7r~ threshold into the physical region (see also [7]). 
What makes this channel particularly apt for an investi- 
gation of the cusp, apart from the enormous statistics col- 
lected by the N A48 / 2 collaboration [3] , is the significantly 
larger branching fraction of tt+tt+tt" compared to 

7r^7r^7r+, such that the perturbation of the decay 
spectrum of the latter is very sizeable. This is in marked 
contrast to two other decays that have been studied sub- 
sequently and that display, in principle, cusp structures 
generated by the same mechanism: Kl ^ Stt^ [2,4,8,9], 
and T] Stt^ [8, 10-14]. In both of these, the weaker cou- 
pling to the charged-pion final state diminishes the cusp 
to a mere 1-2% effect on the decay spectrum. 

In this respect, rj' tjtt^tt^ decays offer a more promis- 
ing candidate for an alternative channel to study the 
cusp.^ Indeed, in the isospin limit BR(77' t^tt+tt") = 
2BR(7^' 777r^7r^), so one expects a sizeable effect on the 
TT^TT^ decay spectrum. In this article, we adapt the non- 
relativistic effective field theory (NREFT) formalism de- 
veloped in [5, 6, 8] to these channels. Compared to i^T ^ Stt 
and T] ^ Stt decays, there is a diflFerent "secondary" final- 
state rescattering channel ttt] to take into account. In ad- 
dition, the formalism has to be slightly amended for par- 
ticles in the final state whose mass difference (M^ — M^^ 
in this case) is not small. This study is very timely with 



regards to the upcoming high-statistics 77' experiments 
at ELSA [16], MAMTC [17,18], WASA-at-COSY [19,20], 
KLOE-at-DA^NE [21,22], or BES-III [23], which are ex- 
pected to increase the data basis on rj^ decays by orders 
of magnitude. 

The outline of this article is as follows. In Sect. 2, we 
present the generalized non-relativistic effective field the- 
ory framework for r]' rjirn decays, define the necessary 
Lagrangians and perform the matching to Dalitz plot as 
well as TTTT and ttt] threshold parameters. Furthermore, we 
comment on the effects of inelastic channels. In Sect. 3, 
we present the result for the decay amplitudes up to two 
loops, including radiative corrections. This section com- 
prises the central result of this study and can be used 
in forthcoming experimental analyses of these decays. We 
turn around the argument and predict the size and specific 
shape of the cusp in Sect. 4. Our findings are summarized 
in Sect. 5. 



2 NREFT for rj' tjittt 

We consider the neutral and charged decay modes 

7?'(P,0^^°(Pl)^°(P2)7?(P3) , 

7/'(P^0 ^^+(PlK"(P2h(P3) . (1) 



^ The cusp in this channel is briefly discussed in the frame- 
work of unitarized chiral perturbation theory in [15]. 



The charged channel only serves as an "auxiliary mode" 
for the cusp analysis. The kinematical variables are defined 
in the usual way Si = (P^/ — Pi)'^, -§1 + 52 + 53 = M^, + 
+ M| + M|, with Mf, i = 1, 2, 3. We will use 

the notation M^± = throughout. 

We shall proceed along the lines of [5] to develop the 
modified non-relativistic Lagrangian framework for our 
calculation. It is set up in such a way that the results 
are manifestly covariant, with the correct analytic struc- 
ture of the decay amplitude fully reproduced in the low- 
energy region. A consistent power counting scheme can be 
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constructed in a very similar manner to the one in [5,8]: 
we introduce the formal non-relativistic parameter e and 
count the pion and rj 3-momenta (in the 77' rest frame) as 
0(e), the kinetic energies Ti = — Mi as O(e^), and the 
masses of the particles involved as 0{1). Let us remark, 
however, that as opposed to the case of TTTr-scattering, 
where the mass difference between the two scattering par- 
ticles amounts to a small contribution in isospin breaking 
at maximum, for 7r77-scattering we count M^j — M-j^ = 0{1). 
We will thus provide a Lagrangian framework that repro- 
duces the low-energy expansion in a standard manner. As 
the momenta of the decay products can be larger than 
in K ^ Stt decays, it is not a priori clear how fast this 
non-relativistic expansion converges for rj' rjTTTr. In this 
sense, the approach is kin to traditional Dalitz plot param- 
eterizations in terms of polynomials: one has to include as 
many terms as necessary to achieve a good description of 
the data, which we expect to provide with the number of 
terms given in the following. 

The loop expansion using this non-relativistic La- 
grangian framework produces a correlated expansion in 
e and nn as well as ttt] threshold parameters, which we 
denote summarily by avrTr and a^r^, respectively, or by a 
in case we refer generically to both. Each two-particle- 
rescattering increases the order of the loop contribution 
by ae [5]. 



2.1 Non-relativistic Lagrangians 

The full non-relativistic Lagrangian can be split up into 
separate parts describing the r]' rjirir tree amplitude 
and the tttt and ttt] final-state interactions: 

The T]' rjiTTT Lagrangian up to O(e^) is given by 
2 

^J2^' (r/'^ {Wr, - Mr^y r/^+^_ + /i.e.) 

- 2W±^+W±^-) + /i.e.) -h . . . , (3) 
where ^, 77, 77' denote non-relativistic pion, 77, and 77' field 



operators, Wa = — Z\, where A is the Laplacian, 

and Gi, Hi are the low-energy couplings in the neutral 
and charged channel, respectively. The ellipsis stands for 
higher orders in the e expansion. 

We start with tttt final-state interactions and consider 
the three channels in (i) (tt^tt^ tt^'tt'^): (00) (00; 00), (x) 
(H — ; 00), (H — ) (H — ; H — ). For the following discussion we 



introduce the notation 



{Vl)^ = {Wn,+iV) , (4) 



forn = a, 6, c, d. Analogous definitions hold for derivatives 
on the T] field, with ^ ^ rj and Vn 
state Lagrangian can be written as [8] 



Vrj- The TTTT final- 



/c=0,± i 



(5) 



where the first part is the free pion propagator and 

A = X^Ci (^^l^l^a^k + h.c}j 

+ x,D,[{^i),{^lr^a^, + H^i{^a)A^br 

-hi^l^l^a^b^h.C.} 
+ XiFi{{^l),^{^lr^a^, + ^l^l{^a),.W 

- 2hi [i^DMr^a^b + ^H{^a)^{^br) 

+ /i.c.}-h... , (6) 



with xoo 1/4, Xx x+_ 1 and hi = sl-^{M^-\-M^-\- 
Ml + ^d)-> where 5- is the physical threshold of the zth 
channel, explicitly /iqo = SM^o, = 3M^ - M^o, = 
2M^. We have omitted P-wave contributions (in the (H — ) 
channel) as they do not contribute in 77' rjirir as long as 
conservation of C-parity is assumed. The ellipsis denotes 
the omission of higher-order terms in the e expansion. 

In the case of ttt] scattering, we consider the chan- 
nels {i) {tjtt'^ rjTT^): (77O) (770;770), (77+) (77+; 77+) (the 
r]7r~ r]7r~ amplitude is identical to t^tt^ r]7r~^ by 
charge conjugation). For the 7r77 Lagrangian we find 



(7) 



where the first term is again the free particle propagator. 
Before giving the explicit form of the interaction piece, we 
define a differential operator 



(8) 



This operator has to be understood as follows: expand 



about the respective thresholds q = (M^ -h M^oY and 



s^^ = (M^ + M,)^ 
s~^Vi^^\r]^a + h.c. = 



^(,t<,t, + e.) 
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then apply the Feynman rules of the theory in momen- 
tum space and resum the result. Note that the differential 
operator does not violate analyticity and unitarity of 
the S-Matrix in the low-energy region, since it is obeyed 
term by term. 

We can now display the interaction piece, 

-AlJ-W^lrj^a^h.c.}^... , (9) 

with hrjO = -h M^o, hrj^ = -h M^, and Arja = 
M^-M^. The ellipsis again stands for higher-order terms 
in e. We do not consider six-particle couplings, since their 
contribution to the amplitude is negligible (see Sect. 2.3). 



2.2 Matching 

We obtain the couplings of by matching to the stan- 
dard Dalitz plot distribution 

\M{x, y)\^ = lA/'p (1 + ay + by^ + dx^ ^ . . .) , (10) 
where A/" is a normalization constant, and 
_ \/3(si - S2) V3{p°2 - Pi) 



2Mn'Qr,' 



Qt]' 



^ {Mr, + 2M^o) [(M^. - Mr,f - S3] ^ 

^ 2Mn>M„oQr,' 

Qr,' = Mr,' -Mr,- 2M^o , 



(11) 



and we have omitted a C- violating term oc x. The particle 
energies in the rj' rest frame are related to the invariants 
Si according to 



2M 



(12) 



For the charged channel we have to replace M^o M^^ 
in (11). Equation (10) is reproduced, up to higher orders 
in X, y, by the following polynomial amplitude: 



M{x,y) =A/'|l + ^^ + i(6-^)7/2 + ^x2 + ...| , 



(13) 



from which one can obtain the matching to the coupling 
constants Gi according to 



(14) 



3d 



In the isospin limit, the charged-channel couplings Hi 
are related to the Gi via Hi = —\/2Gi (in the Condon- 
Shortley phase convention), which we will assume in the 
following. 

The couplings Ci, Di, Ei, Fi are obtained by matching 
to the effective range expansion of tttt and ttt] scattering. 
The partial wave decomposition of the tttt scattering am- 
plitude is conventionally written as [24] 



T\s,t) = 32^ y (2/ + l)tj{s)Pi{z) 



(15) 



^=0 



where tf (s) is the partial wave amplitude of angular mo- 
mentum / and isospin /, Pi{z) are the Legendre polyno- 
mials, and z = cosO is the scattering angle in the center- 
of-mass system. Close to threshold of the pertinent chan- 
nel, one can perform an expansion in the center-of-mass 
momentum q'^^{s) = A(5, M^, M^)/45 (with the standard 



Kallen function \(a,h, c) 
according to 



2{ah- 



- ac - 



he)) 



Re t\ {s) = q^,[al + hUl, + f^ql + 0{ql,) } . (16) 

As only tttt S-waves are considered in the following, we will 
use the slightly simplified notation ao, a2 for the S-wave 
scattering lengths of isospin and 2, and similarly for the 
effective ranges hi and shape parameters fi. 

We use a definition of the TTTy threshold parameters 
strictly analogous to (15), (16), and denote the S- and P- 
wave scattering lengths and the S-wave effective range by 
ao, ai, and h^. These quantities are modified compared to 
the more conventional parameterization used in ttt; scat- 
tering [25]; our ttt; threshold parameters are related to the 
conventional ones ag^, a^^, ^q^^ by 



ai = — a- ' . 



60 = ^":^- (C + 



2M^Mn 



(17) 

We refrain from including the shape parameter of nrj- 
scattering: its contribution is expected to be tiny, and the 
threshold parameters of Trry-scattering are systematically 
smaller than those of TTTr-scattering; c.f. Appendix B. Fur- 
thermore, even the leading threshold parameters of 7r?7- 
scattering are not easy to come by within a sensible error 
range; see [26] and Appendix B. However, the inclusion of 
higher orders in the modified non-relativistic framework is 
straightforward and can be easily performed, should the 
threshold parameters be obtained to greater accuracy. 
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Fig. 1. Non-relativistic two- loop graph involving the six- 
particle vertex. Double lines denote 77^ and 77, the dashed line 
the 7r°. 



We can fix the couplings of the non-relativistic La- 
grangian by matching the amplitude of the effective theory 
to (16). The TTTT and ttt] Lagrangians generate the following 
tree amplitudes: 



Re 7;^^ 
Re T^'^ 



2C,^8D,ql,^32F,qi,- 



(18) 



Z\2 



The matching conditions can be simply read off: for tttt 
scattering we have [5, 8] 



2C00 
2Cx 



N 

3" 

TV 

3" 

TV 



(ao + 2a2)(l-0 , 
(a2 - ao)(^l + I) , 



2C+_ = -(2ao + a2)(l+0 




81^00 



3 ^ 

N . 



262) . 
bo) , 



3 ' 
N 



6 



(19) 



and the matching conditions for the Fi are identical to 
those for the Di with the replacements bi fi/^- Isospin 
breaking in the S-wave scattering lengths has been taken 
into account at leading order in chiral perturbation the- 
ory [27] in (19), ^ = (M^ - Mlo)/Ml, and N = 327r. For 
TTT] scattering we find 

2C^o = 2C^+ = TVao , SDr^o = 8L)^+ = Nbo , 
4Er^o = 4^77+ = 3A/'ai . (20) 



2.3 Six-particle vertices, inelastic channels 

Six-particle interactions can be included in the modified 
non-relativistic framework by straightforward extension 



+ ... , (21) 



where the ellipsis denotes terms with derivative couplings. 
To give a rough estimate of the contribution to the neutral 
channel, we perform a threshold expansion for the diagram 
in Fig. 1, which amounts to the application of the "clas- 
sical" non-relativistic framework [28]. The real part is a 
constant (albeit divergent), which can be absorbed in a 
redefinition of the coupling Go and thus simply amounts 
to a change of the renormalization prescription. The imag- 
inary part of the diagram is given by 



FqGq 



M^oM, 



1/2 



2567r2 (M^ + 2M^o)3/2^''' 



Q2,+o(q3,)^ (22) 




(B) 



Fig. 2. Contributions of the inelastic p^j channel. The double- 
dashed line denotes the the wiggly line a photon; otherwise, 
the line style is as in Fig. 1. 



We can give an estimate for Fq by matching to chiral per- 
turbation theory. At lowest order (assuming isospin sym- 
metry), we find 



Fn = 



^(5M^ 



'tt + } 



(23) 



where F^^ is the pion decay constant. The imaginary part 
of the diagram can be mimicked by allowing Go to have a 
small imaginary part. Numerically, 



ImGo 
Re Go 



10" 



(24) 



which is sufficiently small to assume the coupling con- 
stants of the non-relativistic Lagrangian to be real. 

One may also wonder about effects of the finite width 
of the 77', or the coupling to inelastic channels other than 
r]7T7T. One effect of the finite width is that the mass of 
the T]' obtains an imaginary part; however, in view of 

/Mrj' ~ 2 X 10~^, this is still negligible if we aim at an 
overall accuracy at the percent level at best. The largest 
decay channel of the r]' other than rjirir is p^j [29], which 
could contribute via the diagrams shown in Fig. 2. The 
rescattering vertex in diagram (B) is anomalous, but there 
is no VV rjir^TT^ vertex in the Wess-Zumino-Witten 
Lagrangian [30,31], so we neglect this term. Diagram (A) 
leads to a complex wave-function renormalization factor 
Z, for which we find 



ImZpO^ 



2M2 



: 5 X 10- 



(25) 



This is therefore also an extremely small correction, which 
can furthermore be absorbed in an overall phase. We con- 
clude that, for the purpose of this investigation, we can 
assume real coupling constants and neglect inelastic chan- 
nels. 



3 The decay amplitudes to two loops 

We use the following decomposition of the amplitudes: 



M > + - = A4*^®® 



^tre 



M 
M 



1-loop 

N 

1-loop 
C 



M 
M 



2-loop 

N 

2-loop 
C 



(26) 
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to underline that the tree amphtude of Cr]'^r]7r7r is modi- 
fied by final-state interactions of one, two, etc. loops. 



3.1 Tree amplitudes 

The tree amplitudes are given by 



2 

where Xk = - Mr^, = 1, 2, 3. 



(27) 



3.2 One-loop amplitudes 

For the one-loop amplitudes we find 

A^]v^°°^(5l,52,53) = Bni{s3)J^-{s3) + SAr2 (^3) JqO (^3) 
+ {Bn3{si,S2, Ss)Jrjo{si) + {si ^82)} , 

A^c ''^(^1, 52, ss) = Bci{ss)J+-{ss) + Bc2{ss)Joo{ss) 
+ {Bc3{si,S2, S3) J^+(si) + {si ^ S2)} , (28) 

with the one-loop function 

iqab{Sk) 



Jab{Sk) = 



(29) 



and the polynomials 



Bn3{si,S2, S3) = 2C^o(5i)|go + GiZ'l^ 



+ G2 
+ G3 



17O" 



g^o(gi) 



A 



"^1 



X { Gi + 2G2 ^1"°+ + 2G3 (Xi - Z?". ) } (30) 
for the neutral channel, and 

Bc2(s3) =a(s3)|^G,X^ + G3^g2o(s3)| , 




(A) 




(B) 



Fig. 3. The two distinct topologies appearing at two loops. 
The first one (A) is a simple product of one-loop functions, the 
second one (B) yields a more complicated analytic behavior. 
Here, the double line denotes the 77^ while single lines stand 
generically for any of the particles in the final state 77, tt^, and 



Bc3{si,S2, S3) = 2G^+(si)|fo + HiZ'^X 



H2 

Hs 



A 



S3 - 52 



si 



X + 21^2^1^,+ + 2Hs (Xi - Z^^_) } (31) 

for the charged channel. The following abbreviations have 
been used: 

Q°=rf+P°(+cyd.), Q^ = ^^^%^, 

77' 

^t± = ^[^^^)-Ma. (32) 
Cbc{Sa) = Cbc + 4.Dbcqlc{^a) + l^Fbcqtci^a) , Qx = q+- ' 



Note in particular that there are no higher orders in the 
TTT] mass difference omitted in the P-wave of ttt] scattering; 
with our Lagrangian definition of the P-wave operators, 
the formulae above are exact. 



3.3 Two-loop amplitudes 

The two-loop amplitudes contain diagrams of two distinct 
topologies; see Fig. 3. We find the following: 

>l^^°°P(5i, 52, Ss) = Coo{ss)Bn2{ss)Joo{ss) 

+ Coo{ss)Bni{ss) ^2Cx{ss)Bc2{ss) Joo{ss)J+-iss) 
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Fig. 4. Leading contribution to radiative corrections in rj' 
r]7i^7i^. The wiggly hne denotes the exchange of a Coulomb 
photon between two charged pions (single full lines). 

+ 2a(S3)Sci(S3)J+_(S3) 

+ ^2Cjjo{Si)Bn3{Si,S2, S3) J^o(Si) 

+ 2GoCooCr,oFo{M„o,M^o,Mr„ M^o, Si) 

+ 4HoC^C^oFo{M^, M^, M^, M^o, si) 
+ 4GoC^oPo{Mn, M^o, M^o, M^, si) + (si ^ sa)} 
+ iGoC^oCooFrjiM^, M^o,M^o, M^o, S3) 

+ 8HoCr,+ C:cFr,{Mr„M^,M^,M^,S3) (33) 

in the neutral channel, and 

A^C°°^(S1,S2,S3) = 2C+_(S3)Bci(S3)J^_(S3) 

+ [C, {S3)Bm (S3) + 2C+_ (S3)BC2(S3)] Joo(s3) J+- (S3) 

+ {2C^+(si)Sc3(sl,S2,S3)J^+(si) 

+ moC+-C^+F+{M^, M^, Mrt, M^, si) 
+ 2GoCxC^+F+{M^o,M^o,M^, M^, si) 

+ moC^+F+{Mr„ M„, M^, si) + (si 

+ 8HoCr,+C+-F^{M^, M^, M^, M„, S3) 
+ iGoC^oCa^F^iM^, M^o,M^o, M^o, S3) 



S2)} 



(34) 

in the charged channel. The analytic form of the genuine 
two-loop function F}^{Ma, M^,^ M^^ s^) can be found 
in Appendix A. The representation above is valid to 
(9(a^^e^, a-j^^a-j^ri^'^ ^ ^^r^^^)- smallness of the ttt] 

rescattering effects, this is well justified. 



3.4 Radiative corrections 

In [6], radiative corrections to K ^ Stt decays have been 
discussed in the framework of non-relativistic effective 
field theory. The results obtained there for the decay chan- 
nels Kl Stt^ and Kl tt+tt^tt^ can be adapted im- 
mediately to T]' rjirir decays, so we only quote the final 
formulae and refer to [6] for the derivation. 

In the neutral decay r]' t^tt^tt^, no "external" radia- 
tive corrections exist. The leading electromagnetic con- 
tributions (of (9(a7^7^ loge)) are due to virtual-photon ex- 
change inside a charged-pion loop; see Fig. 4. These be- 
come important very close to threshold, as they modify 
the analytic structure near the cusp by adding a logarith- 
mic singularity to the square-root-like behavior. The dia- 
gram in Fig. 4 can be taken into account by the following 



replacement of the charged-pion one-loop function: 

J+_(S3) ^+-(53) + Jciss) , 

44- (^3) 



Jciss) = --—\og 
oZtt 



Ml 



(35) 



As Tj' r]7T~^7T~ only serves as the "auxiliary chan- 
nel" in the cusp analysis, only radiative corrections of 
O(e^a^) are considered, hence we neglect photon exchange 
inside loops. The external photon corrections, comprising 
virtual-photon exchange as well as real-photon radiation 
up to a maximal photon energy can be subsumed in 
the correction factor i?c(<§3, ^*) that multiplies the decay 
spectrum, which is given in the soft-photon approximation 

by 



dr 

dss 



dss 

2/1 I 



nc{s3,E*) = i + ^l^-^ 

7T I Za 



-a 
3 



log 



(36) 

2^;* 111 

" 3J J ' 



where a = y^l — 4M^/s3. To the accuracy considered 
here, dF^^^/dss is the 77' t^tt+tt" decay spectrum with- 
out photon corrections. See [6] for the more elaborate re- 
sult without the soft-photon approximation. 



4 Prediction of the cusp 

In the previous section, we have given the representation of 
the T]^ T^TTTT decay amplitudes to two-loop order in terms 
of a set of coupling constants G^, Hi, as well as tttt and ttt] 
threshold parameters. In order to extract the tttt scattering 
length combination ao — a2, one would have to fit the latter 
together with the couplings Gi (and possibly even Hi if 
one wishes to relax the assumption on isospin conservation 
in the polynomial terms) to experimental data. Here we 
wish to predict the cusp in rj' tjtt^tt^. For this purpose, 
we first have to fix parameters. 

The TTTT scattering lengths and effective ranges are 
set to the theoretical values ao = 0.220 ± 0.005, a2 = 
-0.0444 ± 0.0010, bo = (0.276 ± 0.006) x M-^, 62 = 
(-0.0803 ±0.0012) X [32], and the shape parameters 
/05 /2 to 0. There is no experimental information on the 
TTT] threshold parameters, and, as discussed in Appendix B, 
theoretical constraints from chiral perturbation theory are 
not very restrictive. We will therefore discuss the un- 
certainties induced by this lack of knowledge by varying 
these parameters in the ranges ao = (0 . . . + 16) x 10~^, 
bo = (0... + 10)x10-3M-2, ai = (-1... + 1)x10-3M-2. 
We find however (in agreement with [33]) that the effects 
of the P-wave di are absolutely negligible (about two or- 
ders of magnitude smaller than the changes induced by 
the variation in the ttt] S-wave shown in the following), 
and therefore set ai = in the sequel. 

There is no very precise information on the 77' tjtttt 
Dalitz plot parameters. We will use the central values 
(without errors) of the most recent determinations done 
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Fig. 5. The decay rate dF/dss divided by phase space. The 
DaUtz plot parameters have been matched to the tree-level am- 
plitude. The dashed line is the tree result, the gray band shows 
the full result under variation of the ttt] threshold parameters. 
The insert focuses on the cusp region around the tt^tt" thresh- 
old. 
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Fig. 6. The decay rate dF/dss divided by phase space. The 
Dalitz plot parameters have been matched to the full ampli- 
tude. We show tree (dashed) and full result (solid), the latter 
again for varying ttt] threshold parameters. Effects of these pa- 
rameters are observed to be significantly reduced compared to 
Fig. 5. The insert magnifies the 7r"^7r~ threshold region. 



by the VES collaboration for the charged channel [34], 
which, adjusting the normalization to the neutral decay, 
read a -0.133, b = -0.116, d = -0.094. We have 
checked using rather different older data on the neutral 
channel (a = —0.116, b = d = [35]) that, while the over- 
all Dalitz plot distribution of course looks very different, 
all statements about the cusp behavior hold in exactly 
the same way. We only discuss normalized decay spectra, 
setting Af = 1. 

The decay spectrum dF/dss for 77' rjir^n^ is calcu- 
lated by 

dr 1 fsti^s) 

d^=512;^i-,3) (37) 



Ml 



S3 



which we normalize by the phase space factor 



4{ss)-Si (S3) 
5127r3M3 



(38) 



In Fig. 5 we display the normalized decay spectrum 
n~^{ss)dr/dss for rj' t^tt^tt^, with the coupling con- 
stants Gi, Hi matched to the VES Dalitz plot parameters 
directly according to (14). The cusp effect is prominently 
visible below the charged-pion threshold and amounts to 
a reduction of the decay spectrum up to around 10%. Fur- 
thermore, the final-state interaction reduces the spectrum 
for large 53 by about 5%, largely due to tttt rescattering. 
The effects generated by varying ttt] threshold parameters 
(as seen by the gray band in Fig. 5) are moderate, but 
particularly small in the cusp region. We remark that the 



sign of the cusp is fixed by (approximate) isospin symme- 
try. Hi = -V2Gi. 

It is obvious from Fig. 5 that the amplitude does not 
reproduce the experimental Dalitz plot parameters any 
more - they are renormalized by the final-state inter- 
actions. We therefore re-adjust the tree- level couplings 
such that the full amplitude squared (10) yields the VES 
parameters. The renormalized tree-level parameters, ac- 
cording to (14), e.g. correspond to ATren = 1.015 (1.021), 

aren = "0.168 (-0.181), 6ren = 3O.IO8 (-0.106), d^en = 

-0.099 (-0.097), for do = 0.016, bo = 0.010 {do = 0, 
bo = 0). The decay spectra resulting from this procedure 
are shown in Fig. 6. We see that the full result follows 
the tree-level spectrum closely except for small deviations 
close to the kinematic limits, and the prominent cusp be- 
low the 7r+7r~ threshold. The uncertainty band due to irrj 
scattering has shrunk to a very narrow line: the effects of 
the third-particle rescattering can be absorbed to a large 
extent in a redefinition of the polynomial part. 

Integrating the spectrum in the region 4M^o < ss < 
4:M^^ we find that the cusp reduces the number of events 
in that region with respect to the tree distribution (or no 
TTTT rescattering) by more than 8%, compared to about 
13% in Tr^TT^TT^ (see e.g. [36]), or less than 2% in 

T] Stt^ [11]. 

In Fig. 7 we zoom in further into the cusp region to 
investigate the expected two-loop cusp above threshold. 
We plot the decay spectrum, with the tree spectrum sub- 
tracted, shifted to at threshold. 



R{ss) = H-\ss) 



■ull 



H-\AMl) 



dss dss 

dPfull dFtree 



ds?, 



ds?, 



(39) 



S3=AMl 
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Fig. 7. Difference between fuh and tree decay rate dP/dss in 
the cusp region, divided by the phase space, and shifted to 
at S3 = 4M^. Note that the scale above the cusp has been 
increased by a factor of 10. 
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Fig. 8. Difference between fuh and tree decay rate dP/dss in 
the cusp region, divided by the phase space. The dashed band 
is without, the fuh band with radiative corrections. The latter 
displays the logarithmic singularity at threshold; compare (35). 



We increase the scale above threshold by a factor of 10. 
Obviously the two-loop cusp (above threshold) is highly 
suppressed compared to the one-loop cusp (below thresh- 
old). One would expect the two-loop cusp effect to be 
smaller than the one generated by one-loop diagrams, as 
it is suppressed by another power in tttt or ttt] thresh- 
old parameters; however, a suppression by more than two 
orders of magnitude may seem surprising at first. The ex- 
planation for this observation can be found resorting to 
the threshold theorem [1,5,6]: in the direct vicinity of the 
7r+7r~ threshold, one can parameterize the amplitude (up 
to three loops, in the absence of photons) as 

Mn = ao^i{ai^ a[q^-{ss)) + (0^2 - 0^2^+- (^3)) 

+ i{as + a'^q+-{ss))+0{a^) , (40) 

where the real parameters ce^, a'- are of 0(a*). This leads 
to 

|A^iv|^ = reg. - 2(ao{a[ + a'^) + a2a[ + ceiceQ 

x^-ql-{ss)^0{a') (41) 

below threshold, where "reg." denotes polynomial (non- 
singular) terms in the vicinity of 53 = 4M^, and 

IMn^ = reg. - 2(ao4 " ^^ic^;)^+-(^3) + 0{a^) (42) 

above threshold. The threshold theorem states that the 
coefficients a[ are proportional to A^c in the appropriate 
kinematics, 

Mc{si = 52,53 = 4M^) oc a; +za2 + a3 + (9(a^) , (43) 
where the factor of proportionality includes the tttt charge- 



exchange scattering length, while obviously 
Mn{si = S2,53 = 4M^) = ao^iai^a2^0{a^) . (44) 

In the isospin limit ai oc —\/2a[_^i and, according to (42), 
the two-loop effects on the cusp would exactly cancel. Nu- 
merically, we find with our set of parameters 

Mc{si = 52, S3 = 4M^) = -1.22 (tree) - 0.076i (1-loop) 

-h 0.0046 (2-loop) , 

Mn{si = S2, 53 = 4M^) = 0.86 (tree) + 0.050i (1-loop) 

- 0.0039 (2-loop) , (45) 

which leads to a suppression of the two-loop cusp by about 
a factor of 250. Incidentally, (41) also allows to estimate 
the effect of three-loop (or 0{a^)) contributions to the 
cusp, which do not vanish in the isospin limit. We find 
that 0{a^) terms ought to reduce the leading 0{a) cusp by 
about 0.5%. Therefore we conclude that, unlike in 
^0^0^+ (jgcayg^ where two-loop contributions are an es- 
sential ingredient to the proper theoretical description of 
the amplitude in the threshold region, the cusp in ry' 
T^TT^TT^ is entirely dominated by the leading 0{a) rescat- 
tering effects. 

Finally, we briefly comment on the effect of the ra- 
diative corrections discussed in Sect. 3.4. Figure 8 shows 
the cusp region similarly to Fig. 7, with and without the 
replacement (35). While the logarithmic singularity at 
threshold is visible, the effect becomes very small away 
from 53 = 4M^. These corrections clearly only matter 
in experimental analyses with very high resolution and 
statistics. 
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5 Summary and conclusion 

In this article, we have generahzed the formahsm of non- 
relativistic effective field theory to describe the analytic 
structure, and in particular the cusp effect, in rj^ rjirn 
decays. We have shown how to construct an effective La- 
grangian that reproduces the next-to-leading threshold 
parameters (effective range, P-wave scattering length) in 
TTT] scattering, and derived the decay amplitudes up to two 
loops, including 0{e^), O(ae^), (9(a^^e^, a^^a^^e^, a^^e^). 
These amplitudes are the central result of our investiga- 
tion, and ought to be employed in future precision studies 
of the rj' rjTTTT Dalitz plot. 

Invoking theoretical information on the coupling con- 
stants involved, we have also predicted the size of the 
cusp effect, and shown that it reduces the decay spec- 
trum below the charged-pion threshold by more than 8%. 
This is a much more sizeable effect than e.g. in rj Stt^ 
decays. Approximate isospin symmetry dictates that the 
cusp of O(a^) above threshold is strongly suppressed, and 
three-loop effects can be estimated to yield a correction 
below 1%. Therefore the threshold singularity in ry' 
rjir^TT^ is determined to very high precision by the lead- 
ing 0{a) rescattering effect. Experimental verifications of 
these predictions at various laboratories [16-23] are ea- 
gerly awaited. 
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A The two-loop function 

The analytic representation of the two-loop formula is 
given as [8] 



V2 



-y2 arctan ■ 



^-y2 

Sk = (Me + Mdf , 

-B =F VB'^ - 4AC 



, y2 = y2{sk) , 



2/1,2 



2A 



B = ql,{sk) -A'^ , C = -qUsk) , 



Sk 



So = M^. + - 2Mrj^ (^M^ + 



2\ 1/2 



4M2 



Sk 



(A.2) 



B TZT] threshold parameters 

Trry scattering has been calculated up to 0{p^) in chiral 
perturbation theory (ChPT) in [25]. We use the form 
quoted in [37] for the expansion near threshold. The 
scattering lengths a^^ and al^ are discussed extensively 
in [26]. 

The S-wave ttt] scattering length to O(p^) is given by^ 



M: 



ao 



(^Ll + U, + ^-Ll-^ + Ll 



1 



-Lr)M2+(L7+^)M, 

+ (15M2 + M^) log ^ + + mA log ^ 



SM^log^ 



■ log — — + A+ + A- '- 



3(M2-M2) "M^ 



A± = {3Mrj ± M^) 



2y/2M^{M^TM^) 

Mrj±M^ 



arctan — , , , 

\^2Mr,{Mr,TM„)J 



(B.l) 



Fk{Ma, Mb, M„ Ma, Sk) = M{2Ah + Bfo) + O(e^) 



with 



(A.1) 



AT: 

fo- 
fi = 

h = 



XyHso,MlMi) 



yi -^V2 - V2 ^ h) , 

{yi{vi - 1) + y2{v2 - 1) - y2{v2 - 1) + /i) , 



4 {vi -^V2 - V2 ^ h) , 
4 
3 



1 + Ql/sk 



-yi arctan - 



1 



1,2 ; 



where we have made use of the Gell-Mann-Okubo relation 



4M|^ 



, in the 0{p ) corrections, ao can be 
related to the / = 2 tttt scattering length 



M2 



167rF2 \ 
L 



+ 



32M 
167r2F2 



2 r 



L\ + Ll 



3 log: 



+ Y W 

+ n log - 

9 fJL 



(B.2) 



^ All 7r?7 threshold parameters are given in the isospin limit. 
Note that isospin breaking in nrj at tree level only affects the 
scattering length at 0{(rnu — TUdY) and is therefore negligible; 
there are no electromagnetic effects at that order. 
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as well diS Af = Fk/F^-1 and Aqmo (4M|^ - 3M^ - 
M^)/{M^ - Ml) [38] in the form of the low-energy theo- 



rem 



ao 



Ml [ 
967rF2 \ 



1 H + -^T^ — TT^A 



M4 ""0 ^ 3 M2 - M2 ""^ 



4 ^ 1 

3 iDTT^i'i 



9M;^-6M2m2 + M^ 
3(M2 - M2) °S 



^ 3M;^-iom2m2-m^ ^ 



+ 7 M, 



(B.3) 



where Aal = +_M2/(167rF^). We find the (modified) 
TTT] effective range bo and P-wave scattering length ai at 

-6o = ^{(xi + L^ + ^-f)(M2 + M2) 

5 + ^)m,m,| 



1 [4 2] 



2M;^(3M2 - 2M^M^ + 3M2) _ A+ 

^ 9(M^-M,)3(M^ + M,) °^M^ 6 



7[M2 + -M^M^ + M2)log^ " 



/X 6(M^ - M,)2 
135M^ - lOSM^M^ + lOTM^M^ - 2UM^M^ 



{54M^{Mr,-M^f) ' 



ai 



-(q + |)«,m.}.3^{4mJ.o.^ 



■log 



(3M^ - 18M^M^ + 7M2) log ^ 
27M^ + lOSM^M^ - 215M*M2 + 214M^M3 

59M2m^ + 2M^M^ - M^l (18M2(M^ -M^)2) 
12M^ ± 3M3m^ -15M^M2(M^ ± M^) - 



2M2(3M^±M^)(M^TM^) 



(B.4) 



As the TTTT scattering amplitude in SU(3) only depends on 
the linear combination 2L\ + L3 [37], it is obvious that 
one cannot formulate low-energy theorems for 60 and ai 
in terms of tttt threshold parameters (or ttK scattering 
lengths [39]). We therefore refrain from recasting these 
results in alternative forms. 

We use two different sets of 0{p^) low-energy con- 
stants for numerical evaluation [40,41]. A major difficulty 





A 


Kei. [4UJ 


Kei. [4ij 


Lirj 1 


lU fto 


7 fK 
( .0 


1 1^ 7 -1- 9Q Q 

10. / zt zo.y 


y.o It 10.0 


9-1-77 


lO^M^^o 





9.9 zb 22.9 


0.4 ± 18.8 




10^ ai 





0.9 ±3.7 


-0.8 ±2.3 








(tttt) 00 


(7r7r)a; 


(7r7r)+- 


10^ ao 




40.9 ±1.7 


-90.0 ±1.7 


70.2 ± 1.7 


lO^M^^o 




38.5 ±2.0 


-118.8 ±2.0 


78.6 ± 2.0 


10^ ai 








19.0 ±0.3 



Table 1. Numerical results for the ttt] threshold parameters. 
The first column "CA" refers to the 0{p^) values, the second 
and third column show the 0{p^) results, evaluated with the 
low-energy constants taken from [40,41]. Errors are obtained 
naively by adding individual uncertainties in quadrature. The 
lower three columns show the tttt threshold parameters for the 
physical channels elastic tt^tt^ scattering (tttt) 00, charge ex- 
change {7V7v)x, and elastic tt^tt" (tttt)-! For details, see text. 



consists in estimating the combined errors, as the vari- 
ous uncertainties for the low-energy constants are strongly 
correlated. The errors quoted in Table 1 are obtained 
by naive error propagation, neglecting any correlations; 
we consider the uncertainties thus obtained significantly 
overestimated. In the case of ao, we also use the low- 
energy theorem (B.3), with Aa^ = 0.0012 ± 0.0010 [32], 
Agmo = 0.196, and Ap = 0.193 ± 0.006 [29]. Table 1 
shows the results for ao, 60, ai thus obtained; in the cases 
of ao and ai, these are consistent with the findings in [26]. 
The conclusion is that chiral symmetry does not make very 
precise predictions for the ttt] threshold parameters. The 
next-to-leading order corrections for ao can be as large as 
the current-algebra value, the magnitude of the effective 
range is very badly constrained, and not even the sign is 
fixed for ai. For numerical evaluation in the main text, we 
decide to vary the ttt] threshold parameters in the ranges 
ao = (0 . . . + 16) X 10-^ 60 = (0 . . . + 10) X IQ-^M-^, 
di = (—1 . . . + 1) X 10~^M~^. We consider these reason- 
able, although not the most conservative limits possible. 
They comprise the values given in [25] and most of the 
parameter ranges discussed in [26]. 

To put these numbers into perspective. Table 1 also 
shows the tttt threshold parameters [32] for the different 
physical channels (corrected for tree-level isospin break- 
ing in the S-wave scattering lengths [27]), in the same 
units. The errors are approximate only and, in case of 
the S- waves, propagated from the dominant / = thresh- 
old parameters. The comparison demonstrates that, even 
within a large uncertainty range, ttt] scattering in general 
is much weaker than tttt scattering and should therefore 
have far less influence on the decay properties of 77' rjirir 
via final-state interactions. We furthermore remark that 
a simple unitarization model [42] (that reproduces the 
ao(980) resonance within reasonable accuracy) indicates 
that the ttt] phase stays below 5° in the range up to 
Vs < Mrj' - M^o ^ 823 MeV, therefore the influence of 
the ao(980) resonance is not yet severe and the effective 
range expansion still applicable. 
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